www.gusucode.com > 小波分析理论与Matlab 7实现 > 小波分析理论与Matlab 7实现/code/Ch13/examp13_22.m
clc; clear; % 产生heavy sine信号和噪声 init = 2055615866; [xref,x] = wnoise(5,11,sqrt(7),init); xref=xref(1:1000); x=x(1:1000); % 画出原始信号和噪声信号 subplot(221), plot(x), title('原始信号') subplot(222), plot(xref), title('噪声信号') % 使用wpdencmp函数进行信号降噪 n = length(x); thr = sqrt(2*log(n*log(n)/log(2))); %计算阈值 xwpd = wpdencmp(x,'s',4,'sym4','sure',thr,1); % 画出小波包降噪信号 subplot(223), plot(xwpd), title('小波包降噪信号') % 比较小波基降噪结果 xwd = wden(x,'rigrsure','s','one',4,'sym4'); % 画出小波降噪信号 subplot(224), plot(xwd), title('小波降噪信号')